The overall time evolution in phase-ordering kinetics 
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I. INTRODUCTION 
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The phenomenology from the time of the quench to the asymptotic behavior in the phase-ordering 
kinetics of a system with conserved order parameter is investigated in the Bray-Humayun model and 
in the Cahn-Hilliard-Cook model. From the comparison of the structure factor in the two models 
the generic pattern of the overall time evolution, based on the sequence "early linear - intermediate 
mean field - late asymptotic regime" is extracted. It is found that the time duration of each of these 
regimes is strongly dependent on the wave vector and on the parameters of the quench, such as the 
amplitude of the initial fluctuations and the final equilibrium temperature. The rich and complex 
' crossover phenomenology arising as these parameters are varied can be accounted for in a simple 

way through the structure of the solution of the Bray-Humayun model. 
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Phase-ordering kinetics is very much an open problem in non equilibrium statistical physics. Its quite simple 
formulation is in striking contrast with the difficulties encountered in its study, due to the strong non linear nature of 
the problem. During the past two decades much attention has been devoted to this field and much progress has been 
made in the understanding of what goes on after a quench inside the unstable region of the phase diagram M. In 
particular, the late stage of the evolution has been studied with great care leading to a detailed knowledge of the scale 
invariant asymptotic time regime. Despite this progress, however, most results recently obtained are still based on 
uncontrolled ad hoc approximations relying on the unproven dynamical scaling hypothesis, while a systematic theory 
allowing for the derivation of late stage properties from first principles is still far to come g] . 
^vq ■ Leaving aside this difficult problem, in this paper we mainly focus on what happens before the scaling regime is 
entered, aiming at a comprehensive description of the overall time evolution in a system with scalar conserved order 
C^l . parameter. Specifically, we will consider the Cahn-Hilliard-Cook (CHC) model or model B in the classification 
of Hohenberg and Halperin Q . The study of the preasymptotic regime HQ has been hitherto somewhat neglected 
because of the absence of scaling and universality. This feature on one hand makes the understanding of the preasymp- 
totics less urgent and on the other makes it more difficult a subject to study than the late stage. As a result, all what 
we have in terms of early stage theory essentially boils down to the linear approximation Q. Nevertheless, it seems 
clear that an understanding of the complex phenomena taking place in the transition from the preasymptotic to the 
asymptotic regime is an important preliminary step in the development of a full theory of phase-ordering kinetics |^,^| . 
In this respect, the effort undertaken in the last few years B|l(l to make an assessment of the relevance of the large- N 
model in the theory of phase-ordering kinetics has turned out to be quite productive. Summarizing briefly, the 
large- N model is the only exactly soluble model presently available for phase-ordering kinetics with conserved order 
parameter. The asymptotic behavior of the exact solution for the structure factor in this model displays multiscaling 
H in place of the standard scaling usually observed in experiments and simulations. Then, the question was raised 
whether multiscaling is the true asymptotic behavior in all cases. The issue was settled by Bray and Humayun (BH) 
by introducing |I^,|l3j a more general model (BH model) which contains the large- TV model as a particular case and 
which displays standard scaling in the asymptotic regime for all finite values of A, where A is the number of order 
parameter components. From this result follows that multiscaling as an asymptotic property holds only if N is strictly 
infinite, while for finite A~ at most it may be observable as a preasymptotic feature. 
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Continuing this line of reasoning, in this paper we carry out a detailed analysis of the BH model arriving to a 
fairly complete picture of the mechanism regulating the structure of the preasymptotic behavior and, what is more 
important, we find that it applies also to the CHC model. The main result then is that what appears to be the generic 
pattern in the overall time evolution is based on the sequence "early linear - intermediate non linear large- N (or mean 
field) - asymptotic fully nonlinear behavior" combined with a strong dependence of the time scales of these regimes 
on the length scale. To be more specific, in wave vector space we identify three regions characterized by greatly 
different durations of the preasymptotic regimes and, furthermore, these durations can be modulated by the choice 
of the parameters of the quench. The three regions are the one around the peak of the structure factor and the two 
to the right and to the left of the peak. The concurrence of all these elements gives rise to a rich variety of behaviors 
in the preasymptotic phenomenology which can be nicely sorted out analytically on the basis of the equation for the 
BH model, but holds true for the full CHC model as well. 

The paper is organized as follows: in Section || we introduce the model and after defining the observables of interest, 
we discuss the linear, mean field and BH approximations. The numerical solution of the BH model is analyzed in 
Section fp . The corresponding analysis of the simulation of the CHC model is presented in Section IV . Finally, 
Section |V| contains concluding remarks. A partial preliminary account of these results was published previously [|14| . 



II. THE MODEL 



We consider a system described by an ./V-component order parameter which evolves according to the CHC equation 
for conserved order parameter 



dt 



ff(x,t) (1) 



where V{<j>) = §0 2 -I- t^O 2 ) 2 , with r < and g > 0, is the mexican hat potential and if is a gaussian white noise, 



2^ 'AN 

with expectations 



< f)(x, t) > =0 
< n a (x, t)rj p (S', f) > = -2T F 5 a pV 2 5(x - x')S(t - t') 



Here Tp is the final temperature of the quench. In the following we shall consider processes with an high temperature 
uncorrelated initial state 

< ${x,0) > = 

<0 Q (x,O)^(f,O)> = AS a0 6(x-x') 1 j 

where A is the size of the initial fluctuations. In the subsequent time development the average order parameter 
remains zero and the observable of interest is the equal time structure factor C(k,t) given by the Fourier transform 
of the generic component of the correlation function < <p a (x,t)<p a (x' ,t) >. 

If the final temperature is below the critical point we know that the structure factor evolves from the initial form 

C(fe,0)=A (4) 

towards the equilibrium form 

lim C(k, t) = M 2 {T F )6(k) + C { T eq) (k) (5) 

where the first term is the Bragg peak associated to ordering, with an equilibrium value M(Tf) for the order parameter 
in the broken symmetry phase, and (k) is the contribution of thermal fluctuations. The problem is to describe 
how this form of the structure factor evolves out of the initial form (^). Parameters of the quench are A and Tp. 

The physical picture born out of experiments and simulations is clear: after a certain initial transient, order is 
established over domains of size L(t) and within domains thermal fluctuations have relaxed to equilibrium. Once 
this arrangement is achieved, memory of the initial condition is lost and the subsequent time evolution (late stage) 
is essentially limited to the coarsening of the ordered domains with thermal fluctuations enslaved in the equilibrium 
form. Phenomenological arguments [Q tell how the size of domains grows in time L(t) = k~^(t) ~ i 1 / 2 , where k m {t) 
is the peak wave vector of the structure factor and z — 3 or z — 4 in the scalar or vectorial case, respectively. At 
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this stage, the structure factor can be split M into an ordering and a fluctuation component which are, to a good 
approximation, decoupled 

C(k,t) =C ord (k,t) + C TF (k,t). (6) 

The ordering component scales with respect to the size of domains 

C ord (k, t) = M 2 {T F )L a {x)F(x) (7) 

with a = d, the space dimensionality of the system, and x = kL(t). The above scaling form is a finite size representa- 
tion of the Bragg peak and reveals that the late stage evolution is a critical phenomenon characterized by the pair of 
exponents (a, z). The final temperature Tp enters only in the prefactor M 2 (Tp) and the exponents are independent of 
both A and Tp. In this sense the late stage is universal and scale invariant, with an attractive fixed point at Tp = 0. 
Most of the recent theoretical work J|,|l5| has been concentrated on the study of the scaling function F(x), obtaining 
successful results through approximations based on the assumption that scaling holds. Much less is known about the 
process before the late stage is entered, when the time evolution is still sensitive to the actual values of the parameters 
(A,Tp), the scaling form (0) does not yet hold and the separation into an ordering and a fluctuation component is a 
much more delicate issue. In the following we address the theoretical study of the preasymptotic behavior. 



A. Linear theory 



The traditional early stage theory is the linear theory obtained by setting g = in Eq. (|l|). It is then straightforward 
to compute the structure factor 

C(k,t) = AC { l) (k,t)+T F C { £(k,t) (8) 

where 

C { Q l) (k,t) =exp[-2fc 2 (fc 2 +r)t] (9) 

and 

C tI & *) = I 1 " ex P [^k 2 (k 2 + r)t] } ■ (10) 

The distinctive features of linear behavior are: 

• exponential growth for k < ko = \f—r. In particular the peak wave vector of the structure factor remains at a 
constant value k m = ko/V^ (except for very short times when Tp > 0, see Appendix) and 

C(k m ,t) ~exp2A4t (11) 



at k 



(V, '"- n ~ ^ 2TpBt 



for Tp = 
for T F > 



(12) 



for k > ko the structure factor relaxes exponentially to 



C(k) - { % Fj( ^2 + r) f °* ^ > I (13) 



This approximation is useful to describe what happens in the very early stage, when the requirement for the linear 
approximation to hold is met. Namely, when the size of the nonlinear term in Eq. ([!]) is small with respect to the 
linear one. For the initial condition this implies A -C Mq = — r/g (small A), where Mo is the equilibrium value of the 
order parameter at zero temperature, at the bottom of the mexican hat potential. Conversely, if A ~ Mq (large A), 
the linear theory does never apply. However, even when it applies at the beginning, at some point this approximation 
is bound to break down due to the exponential growth for k < ko- Furthermore, the form (||) does not describe 
coarsening, since the peak value of the wave vector remains constant in time, and the scaling form (^|) never applies. 
In short, the linear early stage theory does not connect to the late stage theory. This is true also to any finite order in 
g. Therefore, any approximation aiming at a global description of the phase-ordering process must necessarily involve 
some form of infinite resummation of the expansion in powers of g. 
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B. Large- N model 



By analogy with the theory of critical phenomena, one may explore the 1 /N expansion as a method to organize 
infinite partial resummations. To lowest order (large- N model), one obtains from (|lj) the following equation for the 
structure factor 

aC ^' t} = -2k 2 [k 2 + R(t)]C(k,t) + 2k 2 T F (14) 

with 

R(t) = r + gS{t) (15) 

and 

d d k 
(27r)° 



S{t) = I T^C^t). (16) 



This equation can be solved analytically Q yielding the full time evolution from the initial condition (Q) to the final 
equilibrium form (||). However, the time development of C(k,t) obtained in the late stage of this model differs in 
several respects from the expected behavior (||) and (0). First of all, the contribution of the structure factor associated 
to the buildup of the Bragg peak takes the multiscaling form 

C(k, t) ~ [LikmLfl^^Fix) (17) 

where L(t) — i 1 / 4 , k m (t) is the peak wave vector related to L{t) by 

{k m Lf = dlogL + (2 - d) \og(k rn L) (18) 

and x = k/k m (t). Then, the exponent a(x) and the scaling function F(x) do take different forms for the quenches to 
finite temperature with < Tp < T c and for Tp = 0. In the former case 

/ \ , x _ / 2 + (d - 2)i/j(x) x<x* , . 

a(x) = a T (x) = | 2 x > x* ^ ' 

and 

F(x) = ^ (20) 
x 

with x* = y/2 and ip(x) = 1 — (1 — x 2 ) 2 , while in the latter 

a(x) = ao(x) = dip(x) (21) 

and 



F(x) = 1. (22) 

This result shows that the infinite resummation involved in the lowest order of the l/N expansion, although 
producing equilibration and formation of the Bragg peak, it does so in a qualitatively different way from what 
expected in the full CHC theory, since i) the finite size representation ( |l7| ) of the Bragg peak displays multiscaling in 
place of the standard form of scaling (Q) and ii) the fixed point structure regulating the asymptotic scaling properties 
is different as the spectrum of exponents and the scaling function are not anymore independent of the final value of 
the temperature. At Tp = 0, in place of an attractive fixed point, there is an isolated fixed point. Then, there is a 
line of fixed points for < Tp < T c (all sharing the same asymptotic properties) and, finally, there is yet another 
isolated fixed point at T c . It has now been clarified |HJ that the formation of the Bragg peak in the large- iV model is 
associated to a condensation process rather than to the phase-ordering process. Why this ought to induce multiscaling 
is not clear. In any case, this seems to say that the type of resummation involved in the large- N model is not enough 
to produce phase-ordering, neither phase-ordering can be recovered by perturbation theory about the N = oo limit. 
In the following we shall refer to the behavior described above as mean field behavior. 
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In order to expose some features of the solution which are relevant also beyond mean field theory, it is convenient 
to consider the integral form of Eq. (|l^ ) 

C{k, t) = AC ( Q mf) (k, t) + T F C [ ™ f) (k, t) (23) 
which yields the structure factor as the sum of two contributions 

G^ mf \k,t) =exp j-2fc 2 ^'(ft'[fc 2 + i?(0]j (24) 



c -' >ftt) ^i V $S <25) 



coupled together through the self-consistency relation required by the definition of R(t) in Eqs. fll5| ) and ([161). If 
Tp = obviously only Cg ^ \k, t) enters the self-consistency relation and this yields the asymptotic multiscaling 
behavior characterized by oto(x). Instead, if Tp ^ 0, both terms in the right hand side of (^3|) participate in the 
self-consistency relation. Eventually {k,t) prevails and the asymptotic multiscaling behavior (|l7| ) with ar{x) 

is obtained. However, before reaching this regime, there may exist a time interval where the two terms do compete. 
This clearly depends on the relative magnitude of A and Tp. Less obvious is that it should depend also on the length 
scale and that there should be, therefore, a wave vector dependent crossover time t*(k). In order to understand how 
this comes about, a glance at Fig. ([!]) is sufficient. The exponent a(x) is related to the rate of growth of the structure 
factor at some length scale. Roughly, ao(x) and ap{x) indicate how fast, respectively, C^ 1 ^ ik,i) and C^^\k,t) 
grow at the different wave vectors as time goes on. Fig. (Q) shows that ap{x) > ao(x) everywhere, except at x = 1, 
where ao(l) = ot(1) = d. Hence, if A/Tp is sufficiently large for the behavior of CQ m ^(k, t) to be observable at the 
beginning of the quench, this will yield to the true asymptotic behavior of Cj£ ' (k, t) first where the difference 

6a(x) = qt(i) — ao(x) (26) 

is the largest, and then gradually the crossover from Cq" 1 ^ {k, i) to (fc, t) will propagate in time to the length 

scales where Sa(x) is smaller. Therefore, Sa(x) is the key quantity which controls the wave vector dependence of 
t*(k). In the case at hand, the crossover will take place very fast for large values of x and will then propagate towards 
x* from above. When it takes place around x* it occurs also for small values of x, eventually propagating towards 
x = 1 both from the right and the left. Thus the peak of the structure factor is the wave vector region where the 
crossover time is the longest. The rich variety of behaviors generated by Eq. (^) as A and Tp are varied has been 
studied in great detail in Ref. poft . 

C. Bray-Humayun model 

A significant improvement over the large- N model is obtained in the BH model. By combining the gaussian auxiliary 
field approximation of Mazenko pi with the 1 /TV-expansion, a non linear closed equation for the structure factor is 
derived [ITl 



dC(k,t) 
di 



-2k A [k 2 +R(t)]C(k,t) 
k 2 

-2—R(t)D(k, t) + 2k 2 T F (27) 



with 



D &*) = J$fiiJ ^ d C(k-kut)G{ki-k2,t)C{k 2 ,t) (28) 

and R(t) defined by (|l5|) and ([li]). Even if the large- TV model is contained in Eq. (|27]) as a particular case recovered 
in the limit N — > oo, it is not possible to pinpoint the correction made, due to the uncontrolled character of the 
approximation involved in the gaussian auxiliary field method p5|. Neither it is clear how to proceed, at least in 
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principle, in order to improve over Eq. (|27|). In any case, the merit of the BH equation is that in the late stage 
standard scaling is recovered for any finite value of N. This means that, whatever correction is contained in the 
BH model, it is enough to describe phase-ordering rather than the condensation process appearing in the large-TV 
model. Therefore the BH equation may be regarded as the self-contained definition of a basic model for phase-ordering 
kinetics. Even though the model does not allow for a complete explicit solution, the structure of the formal solution 
in conjunction with numerical analysis allows to follow in great detail the overall time development and to uncover 
how the asymptotic scaling behavior is gradually born out of the non universal prcasymptotic behavior. 
Going over to the integral form of Eq. (P^) we find the generalization of Eq. ( |23| ) 

C(k, t) = Atff H) (k, t) + T F C { T B F H) (k, t) + ^C {BH) (k, t) (29) 

with 

C { BH \k,t) =expj-2fc 2 ^ dt'[k 2 + R(t')]j (30) 



C,<f"(M) = -2k' ! dt' C « ^ R{t')D{k,t') (32) 
Jo Cq (k,t'j 



and where again the self-consistency is implemented through the definition of R(t). In general this is a very complicated 
nonlinear integral equation which must be handled numerically. On the other hand, considerations of the type made 
about Eq. ( p3| ) may be extended to the present case and simplifications may be expected if the conditions for one of the 
terms on the right hand side to prevail over the others are realized. First of all, the non universal features dependent 
on the parameters of the quench (A, Tp) are associated to the first two terms, while the asymptotic behavior is due 
to the third one. Taking only this last term into account, BH have been able to extract the scaling behavior which 
obeys the standard form ([?]) with 

{x 2 for x <C 1 

e - cl (x 3 -i) 2 lnJV forx-l (33) 

e- C2X for x > 1 

where c\ and C2 are constants. The exponential decay for large x has been obtained numerically |l7j while the other 
two results are analytical [[12|]1S[| . Before the asymptotics are reached, the behavior of the structure factor is the 
result of the interplay of all the terms in the right hand side of (|2^) . If there exists a time interval where the first 
two terms do actually dominate over the third one, one may anticipate the occurrence of mean field behavior Jl2| , ^3| . 
Furthermore, on the basis of the previous discussion, one may also anticipate that this will be very much sensitive to 
the relative magnitude of A, Tp, 1/N and that there will be a wave vector dependent crossover time. 



D. Multiscaling analysis 

In order to cover all possible behaviors it is quite useful to introduce the multiscaling analysis of the structure 
factor. Both the standard scaling (|?]) and the multiscaling ( p"7| ) behaviors can be written in a unified general form as 

C(k 1 t) = [C 1 (t)] a(x) F(x) (34) 

with x = fc£ 2 (*) and where Ci(t) and C^t) are two growing lengths. As we shall see, the determination of the 
spectrum of exponents a(x) is quite informative. Numerically this is done by considering C(k, t) over subsequent and 
non overlapping time intervals Ti . After taking the logarithm of Eq. ( |34| ) one has 

\ogC{x/C 2 {t),t) = a{x,t) log d{t) -flog F{x) (35) 



G 



and a(x,t) can be computed as the slope of logC(fc,i) versus log£i(f) for the different times in the interval r$, with 
x fixed and k = x/C2(t). In practice this amounts to the measurement of an effective exponent a(x,t) — d Q°agCi{t) 
of the type introduced in critical phenomena for the study of crossovers between competing critical behaviors. Here 
the time t is the parameter driving to criticality and scaling holds when the effective exponent is independent of 
t. Standard scaling corresponds to £i(t) = C-zit) = Lit) = t x / z and a(x) = d, while multiscaling corresponds to 
C\(t) — L(t)[k m (t)L(t)] 2 / d ~ 1 , £2(2) = k^(t) and a(x) genuinely dependent on x. 

^From the structure factor other quantities can be computed numerically giving insight into the preasymptotic 
behavior. In the following, we will consider the position of the peak k m (t) and its height C(k m ,t). Their asymptotic 
behaviors are 

M*)~*~ 1/3 C(k m ,t)~t d / 3 (36) 

for the CHC model, while in the large- N model from ([l7|) and ( |l8| ) follows 

k m (t) = (^) V4 C(km,t) = t^(lo g t)^-^. (37) 

For future reference it is useful to apply the multiscaling analysis to the linear approximation. When Tp = and 
A > from Eq.(j^) one finds 

2k 4 t r 
logC(xk m ,t) = m #s)log L{k m Lf 

log [L(k m Ly/ d X J L 

which gives the effective exponent 



\2/d-l 



log A (38) 



2k 4 t 

namely the spectrum is of the multiscaling form ( pl[ ) with a prefactor linearly growing with time. This explicit 
time dependence shows the absence of scaling. For A = and Tp > the computation of a(x, t) is slightly more 
complicated and is reported in the Appendix. 



E. Numerical solution 

The numerical solution of the CHC equation (|l|) requires a lot of computer resources. This is due to the need of 
reaching long times, in order to study the asymptotic behavior, and of considering large systems, in order to avoid 
finite size effects. The task of analyzing in detail the effective exponent a(x, t) makes things worse because it crucially 
depends on a very precise estimate of the structure factor. This involves the use of very large systems and the average 
over many realizations of the initial fluctuations and of the thermal noise. 

The Bray-Humayun model, since no averaging over realizations is required, has been solved by simple discretization 
of the equation of motion for the correlation function in real space, which involves only laplacians. The size of the 
system has been chosen 1024 2 for d = 2 and 160 3 for d — 3. 

For the scalar case we have resorted, instead, to the Cell Dynamical System (CDS), a cellular automaton which 
reproduces accurately the solution of the CHC equation, while being computationally much more efficient ]Ts| . We 
have considered the version of the CDS algorithm with f(x) = Atanh(a;) and have taken the parameters A — 1.3 and 
D = 0.5, so that M(T F = 0) = ±0.97767. We have considered lattices of size 512 2 for d = 2 and 128 3 for d = 3. 
In the initial state we have considered a uniform distribution of the order parameter between —6/2 and 6/2 so that 
A = 6 2 /12. For each value of the parameters A and Tp the structure factor has been averaged over a suitable number 
of realizations. The values of k m (t), C(k m , t) and a(x, t) have been obtained by interpolating the spherically averaged 
value of C(k, t) with a cubic spline routine. 



III. RESULTS FOR THE BH MODEL 



In this section we present the multiscaling analysis of the structure factor in the BH model. As pointed out above, 
in the asymptotic regime scaling is of the standard type as long as 1/N is not zero, because the non linear contribution 

C^ H \k,t) in Eq. (29) is eventually the largest for all values of x. Here we are mainly interested in the numerical 
study of the preasymptotic regime and in the identification of patterns due to the interplay of the three terms on the 
right hand side of Eq. (|2£ 



7 



A. Tf = 



We start by considering the case of a quench to zero final temperature, for a system with N = 10 3 and A = 10~ 8 . 
The overall evolution of the structure factor is rather complicate and in order to sort out all the features entering it, 
very good quality of data is needed. Therefore, we present the data obtained for d = 2 which are much less noisy 
and are qualitatively the same as those for d = 3, since Tp = 0. With Tp = the competition is reduced to the first 
and the third term in the right hand side of Eq. (^9| ) . Following the reasoning presented in Section [SL this depends 
on the values of A and N and on the rates of growth of the two terms at the different wave vectors. Fig. (|J) is a 
double logarithmic plot of the time evolution of C(k, t) from the very beginning up to some time inside the asymptotic 
regime. Basic features are the prompt appearance of power law tails on both sides of the peak and a more elaborate 
behavior of the peak. In particular, three sharp regimes are identified (early, intermediate and late stage) which are 
clearly separated by almost abrupt changes in the position k m (t) and height C(k m , t) of the peak at times t\ and t% 
(see also Figs, (j^) and (^)). In the early stage, i.e. in the interval (0,ii), there is no coarsening and k m stays constant 
in time with high accuracy (Fig. (||)) as predicted by the linear approximation. That there should be a linear regime 
could have been anticipated from the tiny value of the initial fluctuations and it is further confirmed by the behavior 
of C(k m , t) in Fig. (|J), which obeys well Eq. (|ll|). On the sides of the peak, according to the linear approximation (H), 
there should be an exponential decay in k, both to the right Cq (k, t) ~ e - 2fc4 * an d to the left CQ l \k,t) ~ e _2fc2r *. 
This holds true immediately after the quench, but very soon (i.e. for t -C ti), even if the peak keeps on staying 
in the same place, the exponential decays are replaced by power laws. These are obviously nonlinear effects. The 

( BH) ~ * ( BH) ~~ * 

competition between Cq (fc, t) and C^ (fc, t) leading to the breakdown of the linear approximation follows a non 

( B H) ~* 

trivial pattern with C nl (k,t) taking over first at small and large wave vectors and then closing towards the peak 
both from the right and the left. The propagation of the breakdown of the linear approximation from high to low 
wave vectors has already been discovered and analyzed in detail in the non conserved case pp| |. Here, we find evidence 
that, through a different mechanism, this breakdown follows a yet more rich pattern characterized by the bilateral 
convergence towards the peak. 

For what concerns the nature of the power laws, for small wave vectors we find C(k,t) ~ k 2 implying that the 
asymptotic behavior ( |33"| ) is established almost from the start. For large wave vectors, instead, we find C(k,t) ~ k~ n 
with n ~ 4. Interestingly, even though this non linear feature sets in even before the appearance of the k 2 power law 
for small fc, it is very different from the exponential decay of Eq. ( |33| ) and it lasts quite long before yielding to the 
formation of the true asymptotic behavior. This phenomenon is not related to the interplay of terms in the right hand 
side of (|29|). Rather, it may be related to the formation in the early stage of unstable localized topological defects. 
Such phenomenon is present also in simulations of the CHC equation with n > d and will be studied elsewhere Ej j . In 

(BH) ~ * ( B H) 

any case, the bottom line is that on the tails there is practically no competition between Cq (k, t) and Cl, (k, t), 
as the latter contribution soon prevails, while a long lasting competition takes place around the peak. We will thus 
concentrate on this region. 

Figs. (||) and (|]) show that immediately after the end of the early regime, at the time t\, coarsening begins with 
power law behaviors for k m ~ t 1 ^ and C(k m ,t) ~ t d / z . However, after some time, both these quantities display a 
pronounced anomaly localized around t% whose origin seems to be puzzling. This is especially true considering the 
behavior of k m (t), which obeys the same power law very precisely, with z = 4, before and after the anomaly (see 
Fig. (|J)). For the peak height, while the power law is well obeyed with d/z — 0.53 for t > t%, the behavior in the 
intermediate regime t\ < t < ti cannot be easily fitted to a power law. The insight into what is going on is obtained 
by making a data collapse of the structure factor (Fig. (^)) which gives evidence for the existence of a crossover, in 
the region around the peak, between two different scaling forms. Both the preasymptotic and the asymptotic scaling 
functions are well fitted by the form 

F(x) - e ~ a{x2 - 1)2 (40) 
which describes the large- N and the BH structure factor |l2|,[l8| around the peak. The change in the value of a in 

" I q^j\ — * ( BH) ~ * 

going from one regime to the other means that around the peak Cq (A;, t) prevails over C^ ; (k, t) long enough not 
only to go past the linear region, but also to leave behind the intermediate time interval (£i ,^2) where the dynamics is 
non linear with the self-consistency relation in Eq. ( p9| ) saturated to a good extent by CQ BH \k,t) alone. As a result, 
non linear mean field behavior shows up around the peak. 

The picture of the overall evolution is completed by the behavior of the effective exponent a(x,t). In Fig. (|^), 
referring to times within the early stage, the mechanism of breakdown of the linear approximation described above is 
clearly displayed. For all the curves, except the last one, there are two qualitatively different behaviors for x > x(t) 
and x < x(t), with x(t) decreasing with time. For x < x(t) the effective exponent a(x,t) is well described by the 
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linear approximation (p9\), while for x > x{t) is essentially flat. This corresponds to the onset of the power law to the 
right of the peak in the structure factor in Fig. (^). The retarded appearance of the other power law tail to the left 
of the peak corresponds to the deviation from the linear form of a(x,t) for small x in the last curve of Fi g, (m . The 
timings of the deviations from linear behavior at the different values of x are a consequence of the shape ( J39|) of the 
effective exponent. Breakdown of linear behavior shows up first at large x, where a{x,t) is large and negative and 
later on around x = where a(x, t), although small, is positive. 

Next, the behavior of a[x, t) over the subsequent time history in Fig. (0) shows the formation of very pronounced 
peaks to the right and to the left of x — 1 . These are due to the very fast growth of the power law tails in the early 
stage, as it is evident from the spacing of the curves in Fig. (||). In this time regime the peak is the slowest growing 
mode in the system. It should be noted the multiscaling behavior around x = 1, which is produced by the mean field 
behavior of the peak in the intermediate time regime. Later on the curves for a{x,t) flatten and collapse around the 
constant value a(x) = 2. The disappearance of the time dependence from a(x, t) reflects the onset of standard scaling 
in the late stage. It should be noted, however, that the approach to asymptotic scaling is quite a slow process, since 
a small upward curvature of a(x) persists for a very long time. This curvature was observed in early measurements 
of a(x) in a variety of systems and was correctly interpreted as evidence that the large time behavior, although not 
yet fully asymptotic, was definitely approaching the standard scaling behavior p2^ ]. 

As a test of the above analysis let us consider the expected behavior upon making large (A = 1) the initial 
fluctuations and smaller 1/N(N — 10 6 ). The immediate consequences should then be the total absence of a linear 
regime and a considerable extension of the time interval where Cq (k, t) prevails over C nl (k, t). It turns out that 
in this case it is not crucial to have extremely precise data, so we present results for d = 3. Inspection of the evolution 
of the structure factor (Fig. (Q)) reveals a different morphology from the case of small A and the absence of the linear 
behavior of the peak is quite evident. Furthermore, the behavior of the tails is less dramatic and a power law appears 
only at the largest times for large values of the wave vector. More precisely, looking at the time dependence of k m (t) 
and C(k m , t) in Figs. (^) and (Q) we can recognize an early regime, which lasts up to a time very close to ti, followed 
by a second regime characterized by power laws. The features of the early regime can be explained from the behavior 



of Cq BH> (k, t) for short time. Namely, the exponential in the solution of the large- N model (24) is peaked at 



kl = - ° 2f ' ■ (41) 

For very small t one has R(t') ~ r + gS(0) and therefore fc^i ~ — (r + gS(0))/2. Since in the case considered 
r + gS(Q) > we have k^ < 0. This means that for very short times C(k,t) is monotonically decreasing and the 
maximum is at k = 0. During this short time interval (not shown in our plots) k^ 1 = oo and C(k m , t) = A. As time 

proceeds, S(t) gets smaller and rapidly f£ dt'R(t') becomes negative while fc^ becomes positive and starts growing. 
This time interval corresponds to the decreasing behavior of k^-{t) in Fig. ([3j) for log(t) < 0. Therefore, in the early 

regime the right hand side of Eq. ( p9|) is dominated by the contribution Cq \k, t). 

Later on coarsening begins with power laws both for k m (t) and C(k m ,t). The behavior of the effective exponent 
in Fig. @ shows scaling since a(x, t) becomes time independent. The shape of a(x) obeys quite well the mean field 
multiscaling form (Elf) for x < x* and the standard scaling form a(x) = d = 3 for x > x* . Namely, in the region where 

cto(x) is negative C^ H \k, t) dominates, while in the region where cto(x) is positive, even though Sa(x) = d — cto(x) 

is also positive, due to the large value of A it is Cq (k, t) that prevails, and not just around the peak, but for all 
values x < x* . This is an intermediate scaling regime, characterized by mean held behavior over an extended range of 
wave vectors. Eventually, there must be the crossover to the asymptotic regime with standard scaling over all length 
scales. This occurs for times larger than those considered in our computation. 



B. Tf > 

With a finite final temperature all three terms are present in the right hand side of Eq. (p^). Without going to the 
most general case, for our purposes it is sufficiently illustrative to consider Tp 3> 1/N. In this case, we expect to find, 
after the early stage, a fairly long intermediate regime characterized by the finite temperature mean field behavior 
discussed in Section ||. The eventual asymptotic regime lies beyond the time interval considered in the computation. 
In the opposite limit Tp <ti 1/N, the intermediate regime is expected to be very much similar to the one analyzed 
above with Tp = 0. 

Choosing d = 3, N — 10 6 and Tp — 1CP 3 let us first consider the behavior of the effective exponent a(x,t) 
when A = 0. In this case an early linear regime is expected. This is demonstrated by the set of the first few 



9 



curves in Fig. (|10|), which follow closely the behavior of Eq. (A15) in the Appendix. In particular, the signature of 
linear behavior with finite temperature is given by the relaxation of a(x,t) to zero for x > x* , which corresponds 
to the fast relaxation ( |l3| ) to the time independent value of C(k,t) for k > fco- After the early sta ge, Fig. ( |l0| ) 
displays the onset of the multiscaling regime with a shape of a(x) well represented by ap(x) in Eq. (fl9|). Hence, 

with this choice of the parameters, in the intermediate regime Cj^, H \k, t) dominates the right hand side of Eq. j2S|), 
yielding finite temperature mean field behavior. Changing the size of the initial fluctuations to the large value A = 1 
(Fig. (|ll|)), in the intermediate regime we find a different multiscaling behavior. For x < x* , a(x) is well represented 
by ao(x) = 3ip(x), while for x > x* a(x) ~ 2 suggesting that in this region a(x) = ar(x). Hence, we find a behavior 
of the type described in the discussion of the large- N model in Section O. Namely, the first two terms dominate the 
third one in Eq. (J29J) , and among them Cq (k,t) dominates for x < x* , while C Tp (k,t) dominates for x > x* . 
The limitation of CPU time has not allowed to pursue the time evolution any further. However, on the basis of 
Eq. (^9|) we expect that there will be a second preasymptotic regime with Cy? (k, t) dominating also for x < x* and 
producing the multiscaling behavior with a(x) = ap{x) for all x, before standard scaling with a(x) = 3 is eventually 
established. 



IV. SCALAR CAHN-HILLIARD-COOK MODEL 



In this Section we analyze the overall time evolution of the structure factor for the ordering dynamics of the scalar 
CHC. The main point is that the dependence of the preasymptotic behavior on the parameters of the quench A and 
Tp follows patterns very similar to those just described in the BH model. 

We consider first the zero temperature quench of the d = 3 system with small initial fluctuations A = 8.3 x 10 -8 . 
Apart from quantitative aspects to be specified below, it is evident from the inspection of Fig. ([12]), that the basic 
elements in the evolution of the structure factor are the same as in the analogous plot of Fig. (Q) for the BH model. 
Here too there is an early regime characterized by linear behavior of the peak and rapid growth of power laws on 
the tails. The early regime is abruptly followed by the onset of coarsening which, as in the BH model, displays an 
intermediate regime delimited by an anomaly in the position and the height of the peak, before the eventual late 
stage is entered. It is then interesting to see if, as in the BH model, this intermediate regime can be associated to the 
manifestation of non linear mean field behavior. Actually, that it is so it can be established in an even more clear cut 
way than in the BH model. In fact, the mean field growth exponent is z — 4 while in the scalar system z = 3. These 
are precisely the values we find from k m (t) (Fig. (|l3|)) which gives 1/z ~ 0.250 ± 0.003 in the intermediate regime 
and 1/z = 0.321 ± 0.005 in the late stage. From C(km,t) in Fig. we find consistently d/z = 0.730 ± 0.002 in the 
intermediate regime while no clear power law can be extracted from the data of the late stage. The existence of the 
mean field scaling regime around the peak in the intermediate regime is also demonstrated by the data collapse in 
Fig. (|l^). The fit to the form (^(j) of the scaling function works well for preasymptotic scaling, as in Fig. (|), while 
is inadequate for the asymptotic behavior. Finally, the comparison of the behavior of the effective exponent a(x, t) 
in Fig. (16) and in Fig. (Q) completes the evidence for the close similarity in the overall time evolution in the CHC 
and BH models with the same quench parameters. The difference in the space dimensionality is inessential. There 
are, instead, important quantitative differences, as mentioned above, in various features of the structure factor. Thus 
for k — ► the scalar system has a fc 4 tail while in the BH model one finds the k 2 tail. On the other side, for large k, 
the Porod's tail k~( d+1 *> of the scalar case is replaced by the exponential decay in the BH model and the analytical 
form of the structure factor is manifestly different around the peak in the two models. These differences show that 
certainly the BH model cannot be used for an accurate computation of the scaling function in a scalar system, but 
are immaterial when considering the gross features of the overall time evolution, as we are doing here. Thus, we 
have enough evidence to establish the sequence "early linear -intermediate mean field - fully non linear late stage" 
combined with the dependence of the crossover time on the wave vector, much in the same way as we have found in 
the BH model. We can go one step further and regard this phenomenology, by analogy with the BH model, as the 
manifestation of the competition among different components in the structure factor. If so, the phenomenology should 
respond to the variation of the parameters of the quench with patterns similar to those already observed in the BH 
model. As a matter of fact, taking A = 1 the overall behavior of the structure factor changes from that of Fig. jl^ ) 
to the one in Fig. (17), closely reproducing the change in morphology observed between Fig . (||) and Fig. (||). The 
same observation applies to the behaviors of k m and C(k m ,t). Comparing Figs. Jl3|) and ( |l4|) with Figs. (||) and (^) 
the similarity is very close and it is definitely established upon comparing Fig. (|18[) with Fig. (|9|) for the effective 
exponent a(x,t). Therefore, by making A large a considerable extension of the intermediate regime is observed, as 
in the BH model, and the manifestation of the mean field behavior takes place in the same way, with a(x) ~ ao(x) 
for x < x* and a(x) ~ d — 3 for x > x*. 
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In the scalar case, however, a less rich variety of preasymptotic behaviors can be obtained by varying the parameters 
of the quench, due to much more rigidity in the model, in the sense that N = 1 does not allow to modulate at will 
the strength of the non linear term, which is always large. In particular, because of the upper bound Tp < T c it is 
not possible to realize the conditions for the temperature term to be overwhelmingly larger than the non linear term. 
For this reason the phenomenology observed with Tp > is not significantly different from the one illustrated above 
(see Fig. (0) and Fig. (0)). 



V. CONCLUSIONS 



In this paper we have studied the global evolution of the structure factor from the very beginning of the quench 
down to the fully developed late stage in a system with conserved order parameter. By making a comparative study 
of the BH and CHC model under variation of the parameters of the quench A and Tp we have been able to identify 
generic features in the preasymptotic behavior. Prominent among these are i) the wave vector dependence of the time 
duration of the early and intermediate regimes and ii) the mean field character of the preasymptotic behavior. These 
properties are tightly linked together. The existence of a mean field regime requires, in fact, multiscaling behavior 
which in turn implies a wave vector dependence of the crossover time. On the whole, the complex phenomenology 
of the preasymptotic behavior can be accounted for in a fairly simple manner through the competition of different 
contributions as exemplified in the structure of Eq. (|2^). While for the BH model this is a direct consequence of the 
form of the equation of motion, the applicability of Eq. (|2^) to the evolution of the CHC model is a non trivial new 
result. It should be noted that the preasymptotic phenomenology presented above is, in principle, experimentally 
observable and that it would be interesting to perform a multiscaling analysis of early time experimental data in order 
to look for the crossover from mean field to truly asymptotic behavior. 

On the basis of the results produced, the BH model qualifies for a reference theory of phase-ordering kinetics which 
captures most of the qualitative ingredients entering in the phenomenology of the process. In other words, the model 
can be used as the starting point for theoretical work aiming at an improvement of the approximation in order to 
have a better performance in the computation of asymptotic properties such as the shape of the scaling function. In 
closing, it should be mentioned that the structure of the overall time evolution analyzed in this paper is by no means 
limited to the case of conserved order parameter. The crossover from mean field to asymptotic behavior is expected 
to be a relevant feature of the preasymptotic regime also when the order parameter is not conserved. However, in the 
latter case it is much more delicate to detect than in the former. With conserved order parameter the crossover is 
made evident by the change in the growth exponent from z = 4 to z — 3 and, most of all, by the spectacular change 
from multiscaling to standard scaling. These elements of discrimination are absent in the non conserved case, since 
the crossover does not involve any change in the value z = 2 of the growth exponent and scaling is standard both in 
the mean field and asymptotic regimes. Work in this direction is in progress. 



APPENDIX: 



In this Appendix we compute the form of the effective exponent a(x,t) during the linear regime for finite quench 
temperature and vanishing initial fluctuations. When Tp > and A = the structure factor ( J10| ) takes three possible 
forms depending on the range of the wave vectors considered. 
If 2k 2 (k 2 +r)t < -1 

C(fc,t) = — —^—exp[~2k 2 {k 2 +r)t] . (Al) 

If \2k 2 (k 2 + r)t\ < 1 

C{k,t) = 2k 2 t. (A2) 

If 2k 2 (k 2 + r)t > 1 

C ^ t ) = vr—- ( A3 ) 

K + r 

The boundaries between such ranges are fixed by the conditions 

2k 2 {k 2 +r)t = l and 2k 2 (k 2 + r)t = -1. (A4) 



11 



The first has only one real positive solution 



-r + ^/r 2 + 2/t 



(A5) 



for all times. The second has two real positive solutions 



k 2 (t) 



-r+ yjr 2 -2/t 



k 3 (t) 



-r- yJr 2 -2/t 



(A6) 



but only for t > 2/r 2 . 

Hence, two distinct linear regimes can be identified, depending on whether the time is greater or smaller than 
t = 2/r 2 . For t < t the k > axis is divided into two intervals by ki(t). For k ^> hit) the argument of the 
exponential in Eq. (|l0| ) is large and positive yielding the form ( |A3|) , while it is small for k <C h(t) and the structure 
factor takes the form (|A2|). Notice that hit) is proportional to t / 4 for small t. 

Later on, for t > t , the wave vector axis is divided into four intervals by h(t) > k 2 (t) > k 3 (t) > 0. For k 3> ki(t), 
C(k,t) still has the form (|A3|). In the range kz(t) < k < k2(t), the struc ture factor takes the form (Al). In the two 
other intervals (k < fcs(i) and k2(t) < k < k\(t)) it maintains the form (A2). The values of ka(t) and k 2 (t) coincide 
for t — to, whi le th e y g o asymptotically to zero and to y— r, respectively, for large time 

Using Eqs. and Eq 

and for large x one finds 



5j) the form of a(x,t) is readily computed. For a structure factor obeying Eq. (A3) 



a{x, t) = 



2 log C 2 (t) 
log£i(t) ' 



(A7) 



When Eq. (|A2|) holds one has 



a(x, t) 



logA(t) 



(A8) 



and finally, when Eq. (Al) is valid, neglecting logarithmic corrections and defining ko = yj —r/2 one finds 



(A9) 



Notice that only in the latter case a(x, t) does actually depend on x. 

In order to fully determine a(x,t) the explicit time dependence of the scaling lengths C\{t) = LikmV) 2 !^ 1 and 
£2 (t) = fc" 1 is needed and this requires the determination of the time dependence of the peak position k m during the 
two linear regimes. One can easily see that, for t <C t , k m decreases as t -1 / 4 while for longer times t t it reaches 
the constant value k m — k . Using these expressions we then have for t <C t 



C x {i)=t x ' A 



C 2 (t)=t 



1/4 



(A10) 



while for t ^> to 



AO) = t 



- +1/4 



(d = 2) 



(d = 3) 



(All) 



C 2 {t) = k ~ x = constant. (A12) 

We are finally able to write down the form of the effective exponent in the linear regime with finite final temperature. 
For t < t 

( ,s (2 x < zi(t) = k/h{t) , 

Q ^ = \2 x» Xl (t). ( Al3 ) 

For t > t and d = 2 
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a(x, t) 



while for d 



a iog(fcp«) 

logt 

A log(fcg*) 

* logt 

81og(fc ' 1 ) 



Q lQg(fcgt) 

log(fc- 2 t) 



a(x, i) 



log(fe t) 
log(fcgt) 
log(fc- 2 t) 

log(fc~ t) 



ip(x) 



X < X 2 {t) EE fc/fc 2 (*) 

x 2 (t) < x < x 3 (i) ee fc/fe 3 (i) 
x 3 (t) < x < Xl(t) 
a; » xi(i) 

x <C x 2 (t) 

x 2 (t) < x < x 3 (t) 

x 3 (t) < x < Xi(t) 

X > Xi(t). 



(A14) 



(A15) 



The most relevant feature of these expressions is that a(x, t) vanishes for large x as time grows, as opposed to the 
case with A = where the spectrum predicted by linear theory becomes very large and negative for x > y/2, as time 
grows. 
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FIG. 1. Spectrum of multiscaling exponents in the large- iV model with d=3, Tf = (cto(x)) and Tf > (ar(i)). 



13 



FIG. 2. Structure factor in the BH model with d = 2, N = 10 3 , A = 10~ 8 and T F = 0. Different curves correspond to 
exponentially growing times. The same applies to all other plots of the structure factor. 

FIG. 3. Time evolution of the inverse peak position of the structure factor in the BH model. 

FIG. 4. Time evolution of the height of the structure factor peak in the BH model. 

FIG. 5. Scaled structure factor of Fig. @. 

FIG. 6. Effective exponent a(x,t) for short times in the BH model with d = 2, TV = 10 3 , A = 10~ 8 and T F = 0. The 
different curves are computed for subsequent non overlapping time intervals of duration exponentially growing. The same 
applies to all other plots of the multiscaling spectrum. 

FIG. 7. Effective exponent a(x, t) for long times in the BH model with d = 2, N = 10 3 , A = 1CT 8 and T F = 0. 

FIG. 8. Structure factor in the BH model with d = 3, N = 10 3 , A = 1 and T F = 0. 

FIG. 9. Effective exponent a(x, t) in the BH model with d = 3, N = 10 3 , A = 1 and T F = 0. 

FIG. 10. Effective exponent a(x, t) in the BH model with d = 3, N = 10 6 , A = and T F = 1CT 3 . 

FIG. 11. Effective exponent a(x,t) in the BH model with d = 3, N = 10 6 , A = 1 and T F = 10" 3 . 

FIG. 12. Structure factor in the CHC model with d = 3, A = 8.3 x 10" 8 and T F = 0. 

FIG. 13. Time evolution of the inverse peak position of the structure factor in the CHC model with d = 3. 

FIG. 14. Time evolution of the height of the structure factor peak in the CHC model with d = 3. 

FIG. 15. Scaled structure factor of Fig. (|l|). 

FIG. 16. Effective exponent a(x,t) for long times in the CHC model with d = 3, A = 8.3 x 1CT 8 and T F = 0. 

FIG. 17. Structure factor in the CHC model with d = 3, A = 2.083 and T F = 0. 

FIG. 18. Effective exponent a(x, t) in the CHC model with d = 3, A = 2.083 and T F = 0. 
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